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Abstract 

We analyze lepton flavour violating transitions, leptonic magnetic dipole moments 
(MDMs) and electric dipole moments (EDMs) in a class of models characterized by 
the flavour symmetry A4 x Z3 x U(l) fn, whose choice is motivated by the approximate 
tri-bimaximal mixing observed in neutrino oscillations. We construct the relevant 
low-energy effective Lagrangian where these effects are dominated by dimension six 
operators, suppressed by the scale M of new physics. All the flavour breaking effects 
are universally described by the vacuum expectation values (<&) of a set of spurions. 
We separately analyze both a supersymmetric and a general case. While the observed 
discrepancy 5a^ in the anomalous MDM of the muon suggests M of order of a few 
TeV, several data require M above 10 TeV, in particular the limit on EDM of the 
electron. In the general case also the present limit on BR(fi — > ej) requires M > 10 
TeV, at least. The branching ratios for /i — » ej, r — > ^7 and r — > e~f are all expected 
to be of the same order. In the supersymmetric case the constraint from [i — > ej is 
softened and it can be satisfied by a smaller scale M. In this case both the observed 
<5a^ and the current bound on BR(fj> — > &y) can be satisfied, at the price of a rather 
small value for |($)|, of the order of a few percents, that reflects on a similar value 
for 0i3. 
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1 Introduction 



Aim of this paper is to analyze the predictions of a class of models of lepton masses and 
mixing angles concerning processes which violate the individual lepton number, such as 
H — > e7 and r — > /ry. We will focus on a flavour symmetry that contains the discrete 
group A 4 and is particularly successful in reproducing the lepton mixing angles observed 
in neutrino oscillations [1,2]. Such a mixing pattern is surprisingly close to the so-called tri- 
bimaximal (TB) one [3]. The experimental precision on the atmospheric and solar mixing 
angles has considerably improved in the last few years and now the solar mixing angle is 
known with a la accuracy of about one and a half degrees. Since its first measurement the 
atmospheric mixing angle showed a steady tendency to be very close to maximal, and now 
its experimental error is of about four to five degrees. At present the TB scheme perfectly 
matches the data and it predicts a vanishing angle 9 13 [4]: 
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To a high degree of accuracy the TB mixing pattern can be derived by assuming that 
leptons transform as certain multiplets of the discrete flavour symmetry A4 [5-10], realized 
at a very high energy scale Af. The A 4 group is spontaneously broken at a scale (tp) <C Af 
by a set of scalar multiplets </?, the flavons, whose vacuum expectation values (VEVs) 
are subject to a special alignment. In models that have been explicitly constructed and 
analyzed, this alignment is not arranged by hand, but it is the natural consequence of 
the dynamics of the theory [6-10]. Moreover the TB mixing pattern is expected to be 
modified by corrections of order (tp)/Af <C I. When these corrections are accounted for, 
the mixing angle #13 is no longer vanishing and becomes proportional to (tp)/Af. In these 
models the total lepton number is broken at a large scale Al (possibly related to (</?)) and 
light neutrinos are assumed to be of Majorana type. Depending on the specific way the 
total lepton number is violated (either by higher dimensional operators or via the see-saw 
mechanism), the neutrino mass spectrum can have different properties that can be tested 
in future experiments. 

It would be highly desirable to reveal the characteristic A 4 symmetry pattern in other 
types of observables, not directly related to neutrino properties. Such a possibility becomes 
realistic if there is new physics at a much lower energy scale M, around 1 10 TeV. Indeed 
we have several indications, both from the experimental and from the theory side, that 
this can be the case. For instance, the observed discrepancy in the anomalous magnetic 
moment of the muon, the overwhelming evidence of dark matter, the evolution of the gauge 
coupling constants towards a common high-energy value and the solution of the hierarchy 
problem can all benefit from the existence of new particles around the TeV scale. In this 
paper we assume that such a new scale exists and that the associated degrees of freedom 
do not provide new sources of baryon and/or lepton number violation. 

Therefore, at least four different scales are present in our approach: the lepton number 
breaking scale A^, the scale of flavour breaking Af, the scale introduced by the VEVs of 
the flavon fields (ip) and the new physics scale M. A generic hierarchy among the scales is 



M <C (</?) <C Aj with expected to be comparable to or smaller than Af. 

We will adopt an effective field theory approach, where the dominant physical effects 
of the new particles at low energies can be described by local dimension six operators, 
suppressed by two powers of the new mass scale M and explicitly conserving B and L. We 
can account for the flavour breaking effects by requiring invariance of these operators under 
the flavour symmetry and by encoding the symmetry breaking effects in the flavon fields. 
This approach is strictly related to, and indeed inspired by, that of minimal flavour violation 
(MFV) [11,12], where the flavour group Gf contains, in the lepton sector, SU(3) e cX SU(3)/, 
and where the dimensionless symmetry breaking parameters are the Yukawa couplings 
themselves. While such a choice has the advantage that it can accommodate any pattern 
of lepton masses and mixing angles, it does not provide any clue about the origin of the 
approximate TB pattern observed in the lepton mixing matrix Upmns- F° r this reason, 
here we choose a flavour group that includes the discrete factor v4 4 : Gf = A 4 x Z 3 x U (1)fn, 
a sort of minimal choice to properly account for both, the neutrino and the charged lepton 
mass spectra. Once the choice of the group has been made, we will explicitly build the 
operators of dimension six in the lepton sector. They will contribute to physical effects 
like the anomalous magnetic moments (MDMs) of the charged leptons, their electric dipole 
moments (EDMs) and lepton flavour violating (LFV) transitions like fi — > c-f and r — ► /ry 
0. We will separately treat the general case where no further requirement is enforced and 
the super symmetric case, where additional constraints are present. 

Among the observable quantities that we have analyzed, only the anomalous MDM of 
the muon suggests a relatively small scale M, of the order of few TeV. Apart from this 
indication, in the non-supersymmetric case there are several data that require M above 10 
TeV. First of all the EDM of the electron that represents the strongest bound: M > 80 
TeV. The limits on the r decays r~ — > fi + e~e~ and t~ — > e + fi~fi~ push the scale M 
above approximately 15 TeV. Finally, to satisfy the present limit on BR(fi — > cy), we need 
M > 10 TeV, at least. We also find that the branching ratios for — > c-f, r — > //7 and 
r — > e7 are all of the same order. Given the present limit on BR(fi — > cy) , this implies that 
t — ► jjfy and t — > ej have rates much smaller than the present (and near future) sensitivity. 
In the supersymmetric case the bound from fi — > is softened, due to a cancellation 
occurring in the relevant dipole operator. Depending on \((p)/Af\, whose typical value lies 
in the range 0.001 -j- 0.05, M can be as small as 0.7 -j- 14 TeV. Therefore in a portion 
of the parameter space our model can simultaneously fit the observed discrepancy of the 
anomalous MDM of the muon and respect the current bound on BR(fi —>■ ej). The allowed 
values for ((f)/ Af are rather small, of the order of a few percents at most and this represents 
also the expected range of #13. To tolerate such a low scale M, additional mechanisms are 
required to suppress at least the electron EDM and the rates for r~ — > /i + e~e~ and 
r~ — > e + /i~/i~. 

The paper is organized as follows: in the second section we give a more extended 
overview of this approach, without specifying the flavour symmetry and its breaking pat- 
tern. Then we move to the class of models we are interested in. We will characterize them 

x For a discussion on LFV without a specific flavour symmetry, see for instance [13-16]. For a recent 
review on the subject see ref. [17]. 



in section 3. We will display the corresponding low-energy effective Lagrangian in section 
4. Subsequently we will derive the predictions for leptonic MDMs, EDMs and LFV tran- 
sitions both without (section 5) and with (section 6) the additional constraints implied by 
super symmetry. These sections contain the main results of our paper. They are compared 
to those of MFV in section 7 and summarized in section 8. More technical aspects are 
discussed in appendix A and B. 

2 Low-energy effective Lagrangian for lepton flavour 
violation 

We consider the possibility that the pattern of lepton masses and mixing angles that 
emerged in the last years with the discovery of neutrino oscillations is dictated by some 
flavour symmetry, realized at a very high energy scale Aj [18]. So far all the attempts to 
build a theory of fermion masses based on a flavour symmetry indicate that such a sym- 
metry cannot be exact and that the main features of the fermion spectrum are reproduced 
by specific symmetry breaking patterns. In particular, the closeness of the observed lepton 
mixing to the TB mixing pattern is successfully reproduced in models with spontaneously 
broken flavour symmetries and special vacuum alignment properties. Here we assume that 
the flavour group Gf, which will be specified later on, is spontaneously broken by a set 
of adimensional small parameters ($), |($)| <C 1. These should be interpreted as ratios 
between VEVs of flavon scalar fields (p, transforming non-trivially under Gf, and some 
fundamental scale Af, above which an ultraviolet completion of the theory describes the 
details of the flavour dynamics. The consistency of this picture requires that \{<p)\ C Aj. 
We also assume that the total lepton number (or better, the combination B — L) is violated 
at a large scale A^, so that the light neutrinos get their masses through the low-energy 
operator [19]: 

£ u = ^-(Hh) T Y(Wl) + h.c. (2) 
Al 

where I denotes the three lepton doublets of the standard model (SM) and H is the Higgs 
scalar doublet^. In this scenario, the Yukawa couplings yj (/ = I, u, d) of the SM and the 
matrix Y in L v become functions of (<&): 

y/ = y/«*» , Y = Y{{$)) , (3) 

so that, by treating ($) as spurions transforming under Gf as the corresponding parent 
scalar fields, the whole theory is formally invariant under Gf. Under the assumption that 

2 We adopt two-component spinor notation, so for example e (e c ) denotes the left-handed (right-handed) 
component of the electron field. For instance, in terms of the four-component spinor — (e e c ), the 
bilinears ea u e and e c a v e c correspond to ^)^ v Pttye and ip^PRipe [Pl,r = ^(^TJ 5 )] respectively. We take 
a" = (1,3), ^ = (1,-50, a*" = iK^-u"^), = i^a" -a u a") and g„u = diag(+l, -1, -1, -1), 
where a = (a 1 , a 2 , a 3 ) are the 2x2 Pauli matrices. Here the four-component matrix 7^ is in the chiral 
basis, where the 2x2 blocks along the diagonal vanish, the upper-right block is given by and the 
lower-left block is equal to a^. 



|(3>) I "C 1, the functions and K can be expanded in powers of ($) and only a limited 
number of terms gives a non-negligible contribution §. 

Even when G/ and (<&) are completely specified, as in the case we are going to consider 
in this paper, if and Af are much larger than the electroweak scale, we are left with a 
limited number of tests of the above picture. Of course, there are still several parameters 
that remain to be measured or constrained: neutrinoless double beta decay, the mixing 
angle #13, the type of neutrino mass ordering, the absolute neutrino mass scale and the 
CP-violating phases. However it would be highly desirable to find evidence of the flavour 
symmetry in other types of processes. Such a possibility opens up if there is new physics 
at a much closer energy scale M, around 1 10 TeV. Here we assume that such a new scale 
exists and that the associated degrees of freedom do not provide new sources of baryon 
and/or lepton number violation @. 

Therefore, in the language of effective field theories, the dominant physical effects of 
the new particles at low energies can be described by dimension six operators, suppressed 
by two powers of the new mass scale M and explicitly conserving B and L. If we focus on 
the lepton sector only, the leading terms of the relevant effective Lagrangian are: 

C eff = CKT+e cT H^yi l + C u + i^e cT H^a flu F IMl/ Ml + h.c. + [4-f ermion operators] (4) 

where Ckt stands for the kinetic terms, e is the electric charge and e c the set of SU(2) 
lepton singlets. The complex three by three matrix Ai, with indices in the flavour space, 
is a function of (<&): 

M=M((Q)) (5) 

Invariance under SU(2)x U(l) gauge transformations is guaranteed if is any arbitrary 
combination of B^ Vl the field strength of U(l), and u a W^ vl the non-Abelian field strength 
of SU(2). Here we are interested in the component of this combination along the direction 
of the unbroken U(l) em and we identify with the electromagnetic field strength. We 
can imagine that we derive such an effective Lagrangian from a fundamental theory by 
integrating out two different sets of modes. In a first step we can integrate out the flavon 
fields and the possible degrees of freedom associated with the violation of B — L, thus 
obtaining a complete set of mass terms for all the light particles including neutrinos, 
charged fermions and the quanta with masses around the TeV scale. Subsequently, around 
the scale M — 1 -v- 10 TeV, we integrate out these additional quanta and we generate the 
other operators listed in eq. (jl]). The latter are still invariant under Gf, once we treat the 
symmetry breaking parameters as spurions. 

In this way we can capture all flavour symmetry breaking effects and keep a high degree 
of predictability since the expansion in the small symmetry breaking parameters can be 

3 An exception occurs when the dimension five operator L u arises from the see-saw mechanism. For 
instance, if right-handed neutrinos are present, there can be renormalizablc interactions between them and 
the flavons. By integrating out the right-handed neutrinos, we find that Ax, itself is proportional to the 
mass scale associated with ((p) . In this case the dimension five operator contains a simple pole in (<j>) and 
is no longer local in the symmetry breaking parameters. For operators violating the lepton number L we 
will allow for a generalized power expansion containing poles in ($). 

In our convention, the masses of the new particles are of order gM/Aw, if the underlying theory is 
weakly interacting with a typical coupling constant g. 



truncated after few terms 0. Thereby, the same symmetry breaking parameters that control 
lepton masses and mixing angles also control the flavour pattern of the other operators in 
C e ff. Moreover the effects described by these operators are suppressed by 1/M 2 and not 
by inverse powers of the larger scales A/ and and this opens up the possibility that 
they might be observable in the future. 

In a field basis where the kinetic terms are canonical and the charged lepton mass 
matrix is diagonal, where we will denote vectors and matrices with a hat, the real and 
imaginary parts of the matrix elements jCi u are proportional to the MDMs a, and to the 
EDMs di of charged leptons, respectively [12,16,17,21]: 

v * v * 
di = 2mi—= ReM.u , d; = e—= Im/Au (i = e, u, r) , (6) 

where v ~ 246 GeV is the electroweak breaking scale, defined as a/2 times the VEV of 
the real, electrically neutral, component of H. The off-diagonal elements .My describe the 
amplitudes for the LFV transitions \i — > ej, r — > /ry and r — > ej [12, 13, 16, 17,21-23]: 



BR(k -> Ijj) 12^/2 



n 3 a 



\M lj \ 2 + \M 3l \ 2 ) (7) 



BR(k -> IjViVj) G 3 F m 2 M 4 

where a is the fine structure constant, Gp is the Fermi constant and m 8 is the mass of 
the lepton Z;. Finally the four-fermion operators describe other flavour violating processes 
like r~ — > n + e~e~ and r~ — > e + /i~/i". In this paper our focus will be mainly on the 
processes \x — > e^j and r — > /X7 and we will make only some comments on the four-fermion 
operators. 



3 Models with A4 flavour symmetry 

So far the discussion has been rather general and applies to any theory with a spontaneously 
broken flavour symmetry. Several models have been proposed to produce the TB mixing 
scheme [5-10,24] and among the most economic and simplest ones are those based on the 
discrete group A$ [5-10]. A A is the group of even permutations of four objects, isomorphic 
to the group of discrete rotations in the three-dimensional space that leave invariant a 
regular tetrahedron. It is generated by two elements S and T obeying the relations [25]: 

S 2 = (ST) 3 = T 3 = 1 . (8) 

It has three independent one-dimensional representations, 1, 1' and 1" and one three- 
dimensional representation 3. We present a set of generators S and T for the various 
representations, and the relevant multiplication rules in appendix A. The group A4 has two 
obvious subgroups: G5, which is a reflection subgroup generated by S, and Gt, which is 

5 A non-analytic dependence on the symmetry breaking parameters ($) can be induced by renormaliza- 
tion group running from the high energy scale (ip) down to energies below the electroweak scale. Such a 
dependence can only be determined when all the degrees of freedom lighter than (tp) are known. In MFV 
similar effects have been studied in ref. [20] . 



the group generated by T, isomorphic to Z 3 . These subgroups are of interest for us because 
Gs and Gt are the relevant low-energy symmetries of the neutrino and the charged-lepton 
sectors at leading order, respectively. The TB mixing is then a direct consequence of this 
special symmetry breaking pattern, which is achieved via the vacuum misalignment of 
triplet scalar fields. If <p = (ipi, y9 2 , <Pz) denotes the generic scalar triplet, the VEV 

Moc (1,1,1) (9) 

breaks A4 down to Gs, while 

(^oc (1,0,0) (10) 

breaks A A down to Gt- 

Concerning the flavour group Gt, following ref. [6-8] we will choose 

G f = A 4 xZ 3 x U(1) FN (11) 

The three factors in Gf play different roles. The spontaneous breaking of the first one, v4 4 , 
is directly responsible for the TB mixing. The Z 3 factor is a discrete version of the total 
lepton number and is needed in order to avoid large mixing effects between the flavons that 
give masses to the charged leptons and those giving masses to neutrinos. Finally, U(1)fn 
is responsible for the hierarchy among charged fermion masses [26] . The flavour symmetry 
breaking sector of the model includes the scalar fields <pt, <Ps, £ an d 0. The transformation 
properties of the lepton fields /, e c , //, r c , of the electroweak scalar doublet H and of the 
flavon fields are summarized in table 1. We also assume the following pattern of VEVs for 
the flavon fields: 



(<Pt) 
A/ 

(<Ps) 
A/ 

<0 

A/ 
A/ 



(u, 0, 0) + (cim 2 , c 2 u 2 , c 3 u 2 ) + 0(w 3 ) 
Cb(u, u, u) + 0(u 2 ) 

c a u + 0(u 2 ) (12) 



where (^2,3,0,5 are numbers of order one, while u and t are the small symmetry breaking 
parameters of the theory. All these quantities are in general complex. Notice that we 
have explicitly displayed the sub- leading 0(u 2 ) corrections only for the field (pT- Those for 
the fields (ps and £ are not needed in the present analysis: they have no special pattern 
and are relevant when discussing corrections to neutrino masses and mixing angles. It 
has been shown in ref. [7] how it is possible to achieve this pattern in a natural way, 
as the result of the minimization of the scalar potential of the theory. In that case the 
model was super symmetric, but an analogous pattern can also be obtained in other, non- 
super symmetric, versions. Note that in the supersymmetric case the 0(u 2 ) corrections to 
p> T obey c 2 = c 3 . 

The flavour group of eq. (fTTj) and the set of flavon fields and their VEVs of eq. (fl2"l) can 
be used to obtain neutrino masses with an approximate TB mixing, both in the presence 
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Table 1: The transformation rules of the fields under the symmetries associated with the groups A4, Z3 
and U(l) FiV . 



of right-handed neutrinos from the see-saw mechanism and directly from a set of higher 
dimensional operators [7] (with the only difference that, in the first case, the flavons <ps an d 
£ transform with u 2 under Z3). Therefore the effective Lagrangian, (j3J), and the discussion 
of the related LFV processes, is exactly the same in both cases, at variance with MFV 
where the flavour group in the see-saw case becomes larger [12]. 

Before going into the details of the model, we recall that at the leading order the Yukawa 
couplings in eq. ((3J) are given by: 
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Y = -b/3 26/3 a - b/3 \u , (14) 



where y e , y^, y T , a and 6 are numbers of order one. At this order the mass matrix for the 
charged leptons is diagonal with the relative hierarchy described by the parameter t. We 
will take 

|t| «s 0.05 . (15) 
In the same approximation, the neutrino mass matrix is diagonalized by the transformation: 

v 2 

UTB m vU T B = -r— diag(a + 6, a, -a + 6)w , (16) 
where Utb, up to phases, is the unitary matrix of the TB mixing [3]: 



^/2/3 l/y/3 
Utb = I -1/V6 l/Vs -l/y/2 \ . (17) 

-i/Vq i/V3 +1/V2 



The symmetry breaking parameter u should lie in the range 

0.001 < Id < 0.05 , 



the lower bound coming from the requirement that the Yukawa coupling of the r does 
not exceed 4ir, and the upper bound coming from the requirement that the higher order 



corrections, so far neglected, do not modify too much the leading TB mixing. Indeed, the 
inclusion of higher order corrections modifies all mixing angles by quantities of relative 
order u, especially we should keep the agreement between the predicted and measured 
value of the solar angle within few degrees. The unknown angle 6*13 is expected to be of 
order \u\. Notice that, in most of the allowed range, we have \u\ < \t\. Such a framework 
can also be extended to the quark sector [9,27]. 



4 Low-energy effective Lagrangian with A4 flavour sym- 
metry 



To evaluate the matrix of dipole moments Ai we should analyze the Lagrangian of the 
model: 

£eff = £~KT + C-Y + C-dip + •■• (19) 

Each of the terms on the right-hand side can be thought of as an expansion in powers of the 
flavon fields. Since we have two independent symmetry breaking parameters, u and t, see 
eqs. (112 j) . we consider a double expansion of £ e ff i n powers of u and t. In this expansion we 
keep terms up to the second order in u, i.e. terms quadratic in the fields <ps,T and £. The 
expansion in the parameter t, responsible for the breaking of the Froggatt-Nielsen U(l)irjv 
symmetry, will be stopped at the first non-trivial order, that is by allowing as many powers 
of the 9 field as needed in order to obtain non-vanishing values for all the entries of the 
matrices describing lepton masses and dipole moments 0. Finally, second order corrections 
in u also arise from the sub-leading terms of the VEV (<Pt), eq. (fT2"j) . and will be included 
in our estimates. 



4.1 Kinetic terms 

The expansion of the kinetic terms can be written as: 
The leading order term is given by 

3 3 

C% = tk h^DpU + % J^iKhh^D^ (21) 

8=1 1=1 

where the quantities k and (kfyi can be further expanded in powers of the field 6: 

h = l + k^ + ... , (fcg) 4 = l+*?ff + ... (22) 



6 Concerning the kinetic terms we observe that we can additionally write down operators involving the 
total invariant 6*6 — \9\ 2 . These contribute to the diagonal elements of the kinetic terms. In the kinetic 
terms of the left-handed fields they can be safely neglected, since the leading order correction is of 0(u). 
For the kinetic terms of the right-handed fields, they contribute at the same order as the terms arising 
through a double flavon insertion. 



Note that for our purposes it is enough to expand k up to the zeroth order and (A;g)i up to 
the second order in t. At the first order in u we have only one type of operator contributing 
to the kinetic terms 

u 

(23) 



C { $ r = i-f-{!p T l9 lt D l J) + h.c. 



A 



where k is a coefficient of order one. Here and in what follows, by a parenthesis (• • •) we 
denote an invariant under A±, (• • •)' denotes the 1' singlet and so on. At the second order 
in u we have a richer structure: 



r (2) _ A2) r 



(2) 



KT ~ KTL ~ KTR > 



(24) 



with the labels L and R referring to lepton doublets I and singlets e c , fi c , r c , respectively. 
For lepton doublets we find: 



4?k = *E^ (XJ ^ /} 

i=i f 



(25) 



where X is the list of ^-invariant operators, bilinear in the flavon fields (ps,T and £, 

X = {£ f £,y?T> (^t) 2 ^Wt,^Ws,^^s,^1^ , (26) 

and there are obvious relations among the coefficients hi in order to assure hermiticity. Note 
that the sum in eq. (1251) runs over all bilinears which can couple to form /^-invariants. 
Whether they lead to a non-trivial contribution to the sum depends on the VEVs of the 
flavons. For lepton singlets, we can distinguish a diagonal contribution and a non-diagonal 
one: 

£-KTR = [£-KTR\d + [^KTI^rid (27) 

1 5 

I^ktrU = ij2 E [(KH^)e^D,e c + (k^X^D^ + (h*) i (X i )*°d»D l fl»] (28) 
/ i=i 



= *E 

i=2 
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+ «E 

i=2 
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A3 



-(Xiye^e^D^ + h.c. 



i=2 



■(A' j )"(9t) 2 e=<T"D„T c + ft. 



c. 



(XiYO^^D^T + h.c. 



(29) 



Due to the structure of the flavon VEVs only the term with % = 5 gives a non-vanishing 
contribution in the sum in eq. (129]) . When the flavon fields acquire a VEV according to 
the pattern assumed in eq. (jl2i) . the term Ckt gives rise to non-canonical kinetic terms 
for the lepton fields, of the following form: 



C KT = iKijl^D^ + iKf^D^ 



(30) 



K 




0{u 2 ) 
l + 0(u) 
0{u 2 ) 



0(u 2 ) 
0(u 2 ) 
l + 0(u) 



(31) 



K c 



l + 0(u 2 ,t 2 ) 
0(uH) 
0(u 2 t 2 ) 



0{uH) 
l + 0(u 2 ,t 2 ) 
0{u 2 t) 



0{uH 2 ) 
0(uH) 
l + 0(u 2 ,t 2 ) 



(32) 



As one can see, the corrections to the kinetic terms which render these non-canonical 
are small (at most at order 0(u) and 0(t 2 )). They have to be taken into account when 
calculating the lepton masses and the dipole transitions induced by the matrix M.. How- 
ever, an explicit calculation shows that at the leading order in our expansion parameters 
the charged lepton masses and the elements of Ai are not affected by these non-canonical 
kinetic terms [28] . For this reason we have omitted in K and K c the order one parameters 
that multiply the powers of t and u. 

4.2 Lepton mass terms 

We now move to the operators that give rise to the Yukawa couplings: 

C Y = A + £>» (33) 
There is a part responsible for the charged lepton masses 

A = A (1) +4 2) + - ( 34 ) 



A f A f A f 



+ jte 2 e c H^ (y T ty + (y T ty ' + ^t c h^ (y T i) " + h. c . (35) 



The terms of eqs. (135]) give rise to yi as shown in eq. ffl3|) , if we define 



u — \u e 



y f = y f + y ' f e- 2 ^ (/ = e, n, r) • (36) 



At the next order in u we find: 

4 2) = E ^Vfft «) + £ §^ c ^ f W)' + E W)" + h.c. (37) 



i=2 / i=2 / i=2 f 

where the coefficients y l e , y 1 ^ and y\ are numbers of order one 0. Note that due to the 
special vacuum structure, eq. (TP2~1) . the term with z = 5 in the sum in eq. (1571) gives a 



7 Due to assumed pattern of VEVs, eq. (fT^]) . each bilinear gives rise to a single contribution to the 
lepton masses: out of the several possible ways of combining Xi with I to form the desired singlet, 
only one has a non- vanishing VEV, at order u 2 . Thus there is no ambiguity in the parameters y l f, when 
discussing the lepton mass matrix. 



vanishing contribution. From eqs. (1351) and (1371) . after the breaking of the flavour and the 
electroweak symmetries, we find the following matrix of Yukawa couplings for the charged 
leptons: 



with 



yi 



( y e t 2 u + z ( Pt 2 u 2 
tu 2 



\ 



(2) 2 
Zr 'U Z 



zfh 2 u 2 
(I) 

y^tu + zytu 

(3) 2 
Zr J U Z 



zfh 2 u 2 



\ 



zfhu 2 



(1) 9 

y T u + Z T u J 



(3* 



7 

.( 2 ) 



yfd + y'fC+e-^ + y 2 f + y 3 f e 



+ y 6 f clc b e~ 2 ^ + yjcacte 



■2itp 



y f c 3 + y' f c*e-^ + y^e' 2 ^ + y)c a c\^ 



y f c 2 + y' f cle + y}c' a c b e ^ + y)c a c\^ 



6„* 



-2iip 



-2itp 



for / = e, /i, r. We indicate complex conjugation via * and recognize that for u = \u\e 



iip 



II 



* reads ue~ 21 ^. Note that the off-diagonal elements of yi, proportional to Zf , originate 
either from the sub-leading contributions to the VEV of the (fT multiplet (coefficients 02,3), 
or from a double flavon insertion of the type ^Vs 1 or fas (coefficients c atb ). Note further 
that in the case c 2 = c 3 Zf also equals zf\ 

Similarly, the Lagrangian giving rise to neutrino masses can be expanded as: 

C = 4 1 )+4 2 ) + ... (39) 

with the leading order terms: 

£(1) X a 



At A 



AfA 



(40) 



Due to the non-vanishing Z3 charge of £ and ips the terms equivalent to those in eq. (HOI 



(2) 



with and (p* s are not invariant. Ci leads to eq. (fT4"|) with a = x a c a and b = x b c b . C 
is not relevant for our discussion here, however gives rise to deviations of relative order \u\ 
from TB mixing. 



4.3 Dipole moments 



Coming to the last term of eq. ( fl9i) . the one contributing to the dipole moments, we 
observe that Cdip and Ci have the same structure in flavour space and they only differ by 
the insertion of the electromagnetic field strength a • F = o^F^: 



£d.i 



dip 



Al) £ (2) 
'-'dip ~ '-'dip ~ • ■ ■ 



(41) 



with 



A 1 ) 

dip 



' M 2 
3 



^6 2 e c H^a ■ F (ip T l) + ^O^H^a ■ F ( Vt I)' + ^r c H^a ■ F ((p T l)" (42) 
A f A f A f 



+ ^9 2 e c H^a ■ F ^ T V) + ^d^H^a ■ F (y^/) ' + ^t c H^o ■ F (<^/) " 



ft 



A' 2 



+ h.c. 



and 



'-'dip 



— Y 

M 2 ^ 



i=2 



||#Vtf V ■ F (XJ) + ^Ofrfa ■ F {XJ)' + ^r c H^a ■ F (X,/)" 



/ 



+ h.c. 
(43) 



Since the Yukawa and the dipole couplings involve the same fermion fields and flavons, 
all statements made concerning eqs. (I33|) and (13T1) also hold for eqs. (T4"2"j) and (T4"3"j) . From 
eqs. (|42l and (|4"3l) . after the breaking of the flavour and electroweak symmetries, we find 
the following matrix for the dipole moments: 



with 







( p e t 2 u + 7 ^tV 


7 fW 2 




M = 


V 




7 < 3 V 




= P f + P' f e- 


2iip 








= p fCl + (3' f c\e 




if 


= (3 f c 3 + (3' f c* 2 e 


-«+ + (3}c* a c b e 


^ + /3jc a c 6 e 




= p f c 2 + P' f 4e 







7 < 3) tV 

7 ?V 



\ 



(3 T u + yfPv? J 



(44) 



+ (3y a c b e 

2iip 
2iip 



-2iip 



+ (3}c a cle 



-2iip 



(45) 



for f — e,fi, t. We observe that, in analogy to what occurs in yi, the off-diagonal elements 

(2 3) 

of A4, proportional to 7 j- , originate either from the sub- leading contributions to the 
VEV of the fT multiplet (coefficients 02,3), or from a double flavon insertion of the type 
^fs or £(pg (coefficients c a) &). For C2 = C3 7 j 2 ' 1 and 7 j 3 ' 1 become equal. This is similar to 



.( 2 ) 



.(3) 



in the case of yi. 



4.4 Four-fermion operators 

A complete classification of SU(2) xU(l), Lorentz and Gf invariant four-lepton operators 
goes beyond the scope of this work and here we limit ourselves to some remarks. At variance 
with the operators controlled by the matrix A4, that vanish when either the Higgs doublet 
or the flavon fields are set to zero, four-fermion operators exist that do neither require the 
insertion of the Higgs field, nor the insertion of a flavon field [17,21,23,29]. A first group 
of terms includes 

j^Trrr , ^ff(«) (f,f = e,^r) . (46) 

These operators conserve the individual lepton numbers. When they contain electrons they 
are bound by LEPII data, that typically require a scale M above several TeV. A second 
group consists of: 

^ C ^ C . (47) 

and their conjugates. The second combination in eq. (j4Tj) stands for several independent 
invariants, obtained through different contractions of the involved A± triplets. This other 



set gives rise to interactions that violate the individual lepton numbers L i} with the selection 
rule ALeAL^AL,- = ±2. For instance, one such interaction is given by 



These operators can contribute to LFV decays like r~ — > yU + e~e~, t~ — > e + /i~/i~ and 
their conjugate, whose branching ratios have upper bounds of the order of 10~ 7 [30]. By 
assuming that these operators enter the low-energy Lagrangian with generic order one 
coefficients, through a rough dimensional estimate we find 

1 



G 2 F M 4 



< 1(T 7 (49) 



and we get a lower bound on the scale M of the order of 15 TeV. A more detailed analysis 
of these operators and of the associated phenomenology will be presented elsewhere. 



5 Dipole transitions 

The physical observable quantities we are interested in, MDMs, EDMs and flavour violating 
transitions for leptons, can be computed from the effective Lagrangian of eq. (jl]), whose 
terms have been detailed in the previous section. If we work at the lowest order in the 
symmetry breaking parameter u, the kinetic terms are canonical (zeroth order in u) and y\ 
is diagonal (first order in u). In this approximation also the dipole matrix A4 is diagonal 
and it only contributes to lepton MDMs and EDMs, which we will discuss below. In order 
to find also the size of the flavour violating processes we need to include the sub-leading 
effects originating from insertions of the flavon fields <ps,T and £, shifts of the VEV of 
ifiT and (additional) insertions of 9 in the Lagrangian. As shown above, these generate 
non-canonically normalized K and K c , eq. (1311) and eq. ( 1321) . and render y\ and Ai non- 
diagonal, eq. (|38|) and eq. (f4*4|) . Therefore, we first have to move to the basis where the 
kinetic terms are canonically normalized and the charged lepton mass matrix is diagonal. 
After performing the necessary steps which are outlined in appendix B, we arrive at the 
following leading order result for the dipole matrix that, in this basis, will be denoted by 
M: 

zi 3) yJ T )t 2 u 2 , 

(50) 

zfy r fi T )tu 2 , 



Mee = (3 e t 2 U , 

Me, = ^ (tf%Pe + Vllf - VeZ^P, - 

M er = - 2 {z^yA + f rl f - y e zi 2 ^ T - 

Vr 

M, e = -{y,it ] -^ ] ^)tu 2 , 

Vt 

M Te = -{y Tl ^-z?K)u 2 , 

Vr 



M T » = - (y T j® - 4 3) /i) , 
M. TT = 8 T u . 

The results shown here have been obtained under the assumption that all parameters are 
real, as explained in appendix B, but they remain unchanged in the complex case as far 
as the diagonal elements of M. are concerned. Notice that, in general, all entries of Xi, 
after applying the appropriate transformations for canonical normalization of the kinetic 
terms and diagonalization of the charged lepton mass matrix, are of the same order as the 
corresponding entries of M. before the transformations, see eq. rtHjl . This is expected, if no 
additional cancellation or enhancement takes place. Note further that the structure of the 
leading order result is surprisingly simple, since it does not involve parameters appearing 
in the non-canonically normalized kinetic terms. By having a closer look at the explicit 
form of Xiij one sees that the elements above the diagonal of Xi, Xi efJl , M. eT and Xi^r, 
have a similar structure. Similarly, the elements Xiij below the diagonal, -M Me) Xi Te and 

(2) (3) 

•M r/X , have a common structure. Finally, note that for C2 = C3, which implies Zj = z ) 

(2) (3) 

and 7^ = 7 \ , the elements A4 Te and Ai Tfl become equal. The MDMs and EDMs are 
given at lowest order as: 

a e = 2m e —= Re(3 e t 2 u) , 

V2M 2 V ' 

2m^-^-Re{^tu) , (51) 

= 2m T —= Re(B T u) , 

V2M 2 V ; 



and 



d e = e—= Im(3 e t 2 u) , 

V2M 2 V ; 

df* = e M2 Im0^tu) , (52) 

v 

d T = e—= Im(3 T u) . 

V2M 2 V ; 

As we have anticipated, MDMs and EDMs arise at first order in the u parameter. As a 
consequence their expressions in eqs. (l5Tl) and ( |52l) depend neither on the parameters jjf 

Ik) 

and Zj of the charged lepton mass matrix, nor on the off-diagonal elements of A4, but 
only on /3f, at the leading order. Notice that, for coefficients Bf with absolute values and 
phases of order one, as expected in our model, we have: 

- = 0(2^) , 4 = ( e £) (53) 

We can derive a bound on the scale M, by considering the existing limits on MDMs and 
EDMs and by using eqs. (|53|) as exact equalities to fix the ambiguity of the unknown 
coefficients f3f. We find the results shown in table 2. 

We see from table 2 that, in order to accept values of M in the range 1 -r 10 TeV, we 
should invoke a cancellation in the imaginary part of -M ee , which can be either accidental 



d e < 1.6 x 10 27 e cm 


M > 80 TeV 


d M < 2.8 x 1(T 19 e cm 


M>80 GeV 


Sa e < 3.8 x 1(T 12 


M> 350 GeV 


5a M « 30 x HT 10 


M w 2.7 TeV 



Table 2: Experimental limits on lepton MDMs and EDMs [30-32] and corresponding bounds on the scale 
M, derived from eqs. ([55)) . The data on the r lepton have not been reported since they are much less 
constraining. For the anomalous magnetic moment of the muon, 5a^ stands for the deviation of the 
experimental central value from the SM expectation. 



or due to CP-conservation in the considered sector of the theory. Concerning the flavour 
violating dipole transitions, using eq. (1501 we see that the rate for U — > Zj 7 is dominated 
by the contribution Aiij, since Aiji is suppressed by a relative factor of 0{t) for fi — > e 7 
and r — > /ry and of 0(t 2 ) for r — > e 7 . We get: 

BR(k -> / i7 ) ^ 48vr 3 a 2 



where 



tfy w| (54) 



^ = 4(*r7? ) -«? ) A-) ■ (M) 



The form of w^- shows that for each transition there are two contributions to the rate: 

(2 3) 

one, proportional to 7l , coming from the off-diagonal element (ij) of the original dipole 

(2 3) 

matrix Ai and one, proportional to z\ , coming from the effect of diagonalizing the 
charged lepton mass matrix, see eqs. (J4"4"j) and fl38l) . respectively. Both contribute at the 
same order in it. Since the quantities lOy are all expected to be of order one, we can 
conclude that in the class of models considered here the branching ratios for the three 
transitions fi — > e 7 , r — > /i 7 and r — > e 7 should all be of the same order: 

BRQj, -> e 7 ) « 5i?(r -> /i 7 ) » 5i?(r -> e 7 ) . (56) 

This is a distinctive feature of our class of models, since in most of the other existing models 
there is a substantial difference between the branching ratios [11,12,33,34]. In particular 
it often occurs that BR{ji — > e 7 ) < BR(r — > // 7 ). Given the present experimental bound 
BR(ji — > e 7 ) < 1.2 x 10 _n [35], our result implies that r — > and r — > e 7 have rates 
much below the present and expected future sensitivity [36]. Moreover, from the current 
(future) experimental limit on BR(fi — > e 7 ) [35] and assuming \w^ e \ = 1, we derive the 
following bound on \u/M 2 \: 



u 



BR(p -f e 7 ) < 1.2 x 10~ n (10"") — < 1.2 x 10 -11 (1.1 x 10 - ") GeV" . (57) 

M 2 



Since the parameter \u\ lies in the limited range 0.001 < \u\ < 0.05, see eq. (TT5I) . we find 

\u\ = 0.001 M > 10 (30) TeV (58) 

|«| = 0.05 M > 70 (200) TeV . (59) 

This pushes the scale M considerably above the range we were initially interested in. In 
particular M is shifted above the region of interest for (g — 2)^ and probably for LHC. 

Since Xi Te and M. Tll become equal for c 2 = c 3 , as explained above, also w re and w T ^ in 
eq. (|55|) become equal and therefore also the branching ratio of r — > ey and r — > ^7 will 
be the same, at the leading order. 



6 Supersymmetric case 

As we have seen, the off-diagonal elements of the dipole matrix M. can be traced back 
to two independent sources. They can originate either from >Q*L when the sub-leading 
corrections to the VEV of the fx multiplet (terms proportional to 01,2,3 i n eq. f|T2l ) are 
accounted for or from £ dip , where the relevant double flavon insertions (see eq. (126!) ) are 
considered. In this second case, the only combinations of flavon insertions that can provide 
a non- vanishing contribution are £Vs and its conjugate. In a generic case we expect that 
both these contributions are equally important and contribute at the same order to a given 
off-diagonal dipole transition. There is however a special case where the double flavon 
insertions £Vs and its conjugate are suppressed compared to the sub-leading corrections 
to ifT and an overall depletion in the elements of the matrix Ai below the diagonal takes 
place. This happens, under certain conditions that we are going to specify, when the 
underlying theory is supersymmetric and supersymmetry is softly broken. 

In order to understand the reason why in the supersymmetric case the insertions £'<ps 
and its conjugate are suppressed, it is useful to recall what are the expected sources of 
chirality flip in the lepton sector. In the supersymmetric case we have two Higgs doublets, 
H ud , and the field W of section 2 is replaced by H d . Moreover, in the limit of unbroken 
supersymmetry the lepton mass terms are derived from the superpotential, which is holo- 
morphic in the chiral superfields. At the leading order in the flavour symmetry breaking 
parameters, these terms are of the type: 

Y J d 2 e SUSY T c H d (cp T i)" (60) 

and similarly for the other charged fermions. Therefore the terms in C\ ' containing 
non-holomorphic fields like (p^ and (i = 3, 7) vanish in the supersymmetric limit: 

y' f = , y k f = (k = 3, 7) (61) 

As a matter of fact, in a realistic model supersymmetry is broken, for instance by the VEV 
F of the 2 SUSY component of a chiral superfield. Here we assume that supersymmetry 
breaking originates in a hidden sector of the theory and is transmitted to the observable 
sector via interactions suppressed by the scale Af or by a larger scale, giving rise to soft 



super symmetry breaking terms with a typical mass scale of order msusY- Supersymme- 
try breaking can give rise to non-holomorphic contributions to lepton masses, which are 
however suppressed by powers of msusY / 'A/, as in the following example: 



^2 J d 2 9susYd 2 S usYT c H d (<^ T l)" 9 2 SUSY m S usY (62) 

The ratio msusY / 'A/ is a tiny quantity, which we will neglect in our analysis. 

Beyond fermion mass terms another source of chirality flip is provided by the slepton 
mass terms of left-right type, which involve both an SU(2) doublet and an SU(2) singlet 
slepton. At the leading order in the supersymmetry breaking parameter msusY, such soft 
breaking masses do not contain any insertion of anti-holomorphic fields like tptp s and £t. 
Indeed insertions of holomorphic supermultiplets only require one power of msusY as in 
the following example 



A~ J d 2 S usYT c H d {ip T l)'' 6 2 SUSY m S usY 



(63) 



On the contrary the insertion of an anti-holomorphic supermultiplet requires two powers 
of m S usY- 

^2 J d 2 0susYd 2 9susYT c H d (^ T l)" 2 SUSY 9 SUSY m 2 SUSY . (64) 

In our approximation we will also neglect all contributions of the type in eq. (j64j) . Further- 
more, we assume that in the underlying fundamental theory the only sources of chirality 
flip are either fermion masses or sfermion masses of left-right type. Both of them, up to 
the order u 2 , are described by the insertion of or V^r i n the relevant operators. This 
is our definition of the "supersymmetric case". Therefore in the supersymmetric case we 
should use the restrictions shown in eqs. ( l6Ti) . The fact that any chirality flip up to the 
order u 2 necessarily requires the insertion of or tp\ also applies to the operators of £ d i P , 
and we have (analogously to eq. (l6Tj) ): 

P' f = , P k f = (k = 3, 7) . (65) 

By applying eqs. fIBTl) and (1651) we find the following matrices yi and M. in the supersym- 
metric case: 

y e t 2 u + (y e ci + y 2 ) t 2 u 2 y e c 3 t 2 u 2 y e c 2 t 2 u 2 

Vi=\ C2 tu 2 y^tu + (y^ ci + y 2 ) tu 2 c 3 tu 2 \ (66) 

y T c 3 u 2 y T c 2 u 2 y T u + (y T c x + y 2 ) u 2 



and 



p e t 2 u + {(3 e c 1 +p 2 ) t 2 u 2 (3 e c 3 t 2 u 2 (3 e c 2 t 2 u 2 

M = | (3^ c 2 tu 2 (3^ tu + (^ ci + (3D tu 2 c 3 tu 2 

(3 T c 3 u 2 (3 T c 2 u 2 (3 T u + (f3 T ci + (3D u 1 



(67) 



The procedure to arrive at M. is exactly the same as before. Thereby, note that the form of 
the kinetic terms K and K c is the same as before, since supersymmetry puts no restriction 
on these terms. The result for the entries of Ai is therfj: 



M 
M 

M 

M 

M 
M 

M 
M 
M 



(if I 



P e t 2 u , 

C 2 + C 3 . n \ .2 2 
{Pe Vfi -yePfijt U , 

c 2 + c 3 



{PeVr ~ y e (3r)t 2 U 2 

(ylPp-Plv^tu 3 



fie 



fifi 



flT 



Tfl 



y 

y, 

Pptu 

c 2 + c 3 
yr 

Vr 

C ^{y 2 T p T -p 2 T y T )u 3 
y T 

f3 T U . 



(68) 



(Pi* yr ~ Pr Vp) t U 2 



Again, the leading order results are not affected by the non-canonically normalized kinetic 
terms and therefore are rather simple. Also here, the elements above the diagonal A4y have 
a common structure, i.e. they are determined by the sum of c 2 and C3 and the parameters 
Df and /3y where indices / and /' are determined by i as well as j. Similarly, for the 
elements M.^ below the diagonal we find that they are either proportional to c 2 or c 3 and 
the index / of all parameters is solely given by the index i of M.ij. In this case the equality 
of Ai Te and M. Til for c 2 = c 3 also holds and is more obvious than in the general case. 
The EDMs and MDMs are similar to those of the general non-supersymmetric case: we 
find the same degree of suppression in the t and u parameters for the physical quantities. 
The difference between the general and the supersymmetric approaches becomes manifest 
when the sub-leading corrections are taken into account, i.e. only in the study of the LFV 
processes. The main difference is the suppression of the elements below the diagonal which 
are not of the order as expected from eq. (JBTj) . but suppressed by an additional power of 
u. We get 

BR(li Ifi) 48vr 3 a 



BR(li — > IjViVj) 



G 2 F M 4 



m 



+ 



■ , ( 2 ) |2 



(69) 



where 



^ fie 




-Ply,) 


w w = 

fie 


c 2 + c 3 

y^ye 


(Pey,, 


- ye Pfi) 


(70) 


^Tfl 


y; 


-PrVr) 


w w = 

Tfl 


c 2 + c 3 

yryfi 


(Pfiyr 


- Pr Vp) 


(71) 


re 


= %{ylPr 


-Plyr) 


w m = 

re 


c 2 + c 3 

y T y e 


(p e yr 


- ye Pt) 


(72) 



8 Note that if we take the formulae given in eqs. (|50|) and use eqs. (|6T|) and (|65p we also arrive at this 
form of Aiij. Thereby, the elements Ai^e, -M Te and M. Tp _ below the diagonal turn out to be zero, i.e. we 
can also deduce the suppression of these elements, however not their exact form. 
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Figure 1: The branching ratio of fi — > ej as a function of |tt|, eq. (75). The deviation of the anomalous 
magnetic moment of the muon from the SM expectation is kept fixed to its experimental range [32] . The 

( 1 21 

unknown coefficients w^e are equal to 1 (darker region) or are random complex numbers with absolute 
values between zero and two (lighter region). The continuous (dashed) horizontal line corresponds to the 
present (future expected) experimental bound on BR{^jl — > ey) [35]. 



Notice that now the first contribution on the right-hand side of eq. (|69|) is suppressed by 
a factor of u compared to the non-supersymmetric case. In most of the allowed range of it, 
the branching ratios of /i — > ej and r — ► /ry are similar and larger than the branching ratio 
of r — > ej. Assuming | = 1, the present (future) experimental limit on BR(fi — > ej) 

implies the following bounds 



| it | 
I it I 



0.001 
0.05 



M > 0.7 (2) TeV 
M > 14 (48) TeV 



(73) 
(74) 



We see that at variance with the non-supersymmetric case there is a range of permitted 
values of the parameter \u\ for which the scale M can be sufficiently small to allow an 
explanation of the observed discrepancy in a M , without conflicting with the present bound 
on BR(fi — > erf). We can eliminate the dependence on the unknown scale M by combining 
1531) and ( |69l) . For fi — >• we get@: 



eqs. 



BR(n -> er/) 
BR(fi — > ev^pe 



12?r 3 a . , 2 

5R W 



I 12| |4 , m e 



7,(2) 1 2 



(75) 



(12) 

where w^ e ' are unknown, order one coefficients. We plot BR(fi — > ej) versus \u\ in fig. 
1, where the coefficients w^e are kept fixed to 1 (darker region) or are random complex 
numbers with absolute values between zero and two (lighter region). The deviation of the 
anomalous magnetic moment of the muon from the SM prediction is in the interval of the 



9 Relations among BR(p — > ey), Sa^ and neutrino mass parameters have also been discussed in [37]. 



experimentally allowed values, about three sigma away from zero. Even if the ignorance 

~ (1 2) 

about the coefficients does not allow us to derive a sharp limit on \u\, we see that the 
present limit on BR(fi — ■> ej) disfavors values of \u\ larger than few percents. We recall 
that in this model the magnitudes of \u\ and #13 are comparable. 



7 Comparison with Minimal Flavour Violation 

It is instructive to compare the previous results with those of MFV [11, 12]. If we restrict 
ourselves to the case where right-handed neutrinos do not affect the flavour properties of 
the theory, the flavour group of MFV is Gf = SU(3) e c x SU(3)i x where we have only 
displayed the part relevant for the lepton sector. Electroweak singlets e c and doublets I 
transform as (3,1) and (1,3), respectively. The flavon fields or, better, their VEVs are 
the Yukawa couplings y% and Y themselves, transforming as (3,3) and (1,6), respectively. 
By going to a basis where the charged leptons are diagonal, yi and Y can be expressed in 
terms of lepton masses and mixing angles (we keep using a notation where a hat denotes 
matrices in this particular basis): 

y l = ^ m f^ , Y = ^U*mf^ , (76) 

where U is the lepton mixing matrix. The diagonal elements of the matrix M. evaluated in 
MFV are analogous to those of the previous class of models and similar bounds on the scale 
M are derived from the existing data on MDMs and EDMs. The off-diagonal elements are 
given by: 

Mij = my j Yh 

= V2P^-^[Am 2 sol U l2 U; 2 ±Aml tm U l3 U; 3 } (77) 

where f3 is an overall coefficient of order one and the plus (minus) sign refers to the case 
of normal (inverted) hierarchy. We see that, due to the presence of the ratio A\jv 2 the 
overall scale of these matrix elements is much less constrained than in the previous case. 
This is due to the fact that MFV does not restrict the overall strength of the coupling 
constants Y, apart from the requirement that they remain in the perturbative regime. 
Very small or relatively large (but smaller than one) Y can be accommodated by adjusting 
the scale Al- On the contrary this is not allowed in the case previously discussed where 
the size of the symmetry breaking effects is restricted to the small window ffT8l) and the 
scale Al is determined within a factor of about fifty. The conclusion is that in MFV the 
non-observation of — > ^7 could be justified by choosing a small A^, while a positive 
signal in /1 — > with a branching ratio in the range 1.2 x 10~ n 4- 10~ 13 could also be 
fitted by an appropriate A^, apart from a small region of the #13 angle, around 6*13 ~ 0.02 
where a cancellation can take place. 

The dependence on the scale Al is eliminated by considering ratios of branching ratios. 
For instance: 



BR(fj, — > ej) BR(t — > hVtV^ 
BR(n -> ev^Ve) BR(r -> /ry) 



2Am? 



2 



<1 , (78 



where we took the TB ansatz, see eq. (p^), to fix 9 12 and 23 . We see that BR(fi — > ej) < 
BR(t —>■ fij) always in MFV. Moreover, for # 13 above approximately 0.07, BR({i — ► ej) < 
1.2 x 10~ n implies BR(r — > 717) < 10~ 9 . For #13 below 0.07, apart possibly from a small 
region around #13 ~ 0.02, both the transitions fi — > cy and r — > fi^y might be above the 
sensitivity of the future experiments. 

We also observe that in MFV the only difference between the general case and the 
supersymmetric one is the presence of two doublets in the low-energy Lagrangian. In MFV 
a chirality flip in leptonic operators necessarily requires the insertion of the matrix yi, both 
in the general and in the supersymmetric case and, apart from the possibility of tan/3 
enhanced contributions, similar predictions for the LFV processes are expected in the two 
cases. 

8 Conclusion 

Flavour mixing within lepton families is surprisingly close to the TB one [3] , with a striking 
agreement for the solar and the atmospheric mixing angles. The closeness of the lepton 
mixing to the TB one is still waiting for a full experimental confirmation, which could be 
provided by the evidence for a relatively small angle #13 , of the order of few degrees. In the 
recent years several models of lepton masses and mixing angles predicting an approximate 
TB mixing have been proposed [5-10,24]. A particularly simple model is the one based on 
the flavour symmetry A4 x Z% x U(1)fn [6-8], which represents a sort of minimal choice 
to properly account for both, the neutrino and the charged lepton sector. The TB mixing 
originates from the misalignment of the scalar multiplets tp that break v4 4 giving masses to 
neutrinos and to charged leptons. In this paper we have analyzed a set of new low-energy 
predictions of models invariant under the flavour symmetry x Z 3 x U(1)fn- These 
predictions involve leptonic MDMs, EDMs and LFV transitions like fi — > ej, r — > ^7 
and r — > ej. We have constructed the effective low-energy Lagrangian that describes 
these observables. Such an effective Lagrangian is invariant under the flavour symmetry, 
and all flavour breaking effects are encoded in the dependence on the spurions $ = y?/A/ 
(Af being the cut-off of the theory). These also control lepton masses and mixing angles. 
The dominant operators are obtained by expanding the Lagrangian in powers of $ and by 
keeping the first few terms in the expansion. The leading contributions have dimension 
six and are suppressed by two powers of a new scale M. This scale can be considerably 
smaller than the other scales of the problem, like the cut-off Af, the VEV of the flavons (if) 
and the scale of lepton number violation A^. We treat M as an independent parameter, to 
be constrained by the experiments. Apart from the scale M, all the relevant information 
needed to predict MDMs, EDMs and LFV transitions is contained in a dimensionless 
matrix A4, whose elements can be computed up to unknown order-one coefficients from 
our Lagrangian. The strongest bound comes from the EDM of the electron [31] : M > 
80 TeV. A lower value for M can be tolerated in the presence of a cancellation in the 
imaginary part of Xi ee , perhaps related to CP-conservation in the considered sector of 
the theory. This problem is also present in MFV [11, 12], where one simply assumes that 
M. ee is real. Another stringent bound can be derived from the analysis of the r decays 



t~ — > fi + e~e~ and r~ — > e + fi~fi~ [30] that, barring cancellations in the coefficients of 
the relevant operators, push the scale M above approximately 15 TeV. Coming to LFV 
dipole transitions, we have found that in the general case the branching ratios for /j, — > cy, 
t — > ji^f and r — > e7 are all expected to be of the same order, at variance with MFV. 
Given the present limit on BR(fi — > ej) [35], this implies that r — > 7x7 and r — > e7 
have rates much smaller than the present (and near future) sensitivity [36]. The absolute 
values of these branching ratios depend on the spurion VEVs (<&) that in our class of 
models are determined by a parameter \u\ in the range 0.001 < \u\ < 0.05. In the general 
case, for BR(fi -> e7) < 1.2 x lO" 11 (10~ 13 ) we get M > 10 (30) TeV if |u| = 0.001 and 
M > 70 (200) TeV for |u| = 0.05. The anomalous MDM of the muon a M and its deviation 
from the SM expectation [32] provide the indication for a lower scale M, of the order of few 
TeV, which would also be of great interest for LHC. In order to reconcile this possibility 
with the results derived from the LFV dipole transitions, we have reconsidered the matrix 
Xi in a supersymmetric context, where additional constraints have to be applied. The 
operators describing \i — > ej, r — > 7x7 and r — > e7 flip the lepton chirality. By assuming 
that in a supersymmetric theory the only sources of chirality flips are the fermion masses 
and the sfermion mass terms of left-right type, we find that a cancellation takes place in 
the elements of Ai below the diagonal. As a result the limits on the scale M become less 
severe. For BR(fi -> cy) < 1.2 x 10~ n (10~ 13 ) we get M > 0.7 (2) TeV if |u| = 0.001 and 
M > 14 (48) TeV for \u\ = 0.05. At variance with the non-supersymmetric case there is a 
range of values of the parameter \u\ for which the scale M can be sufficiently small to allow 
for an explanation of the observed discrepancy in a M , without conflicting with the present 
bound on ji — > ej. Since in our framework 6*13 is comparable to \u\, the present limit on 
BR(/i — > ej) together with the existing discrepancy in a M point to a rather small value for 
#13, of the order of few percents in radians, close to but probably just below the sensitivity 
expected in future experiments at reactors or with high intensity neutrino beams [2,38]. 
It remains to be seen if the bounds on the scale M derived from four-lepton operators can 
be evaded in the supersymmetric case. 
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A The group A4 

The group A4 is generated by two elements S and T obeying the relations [25]: 

S 2 = (ST) 3 = T 3 = 1 . (79) 

It has three independent one-dimensional representations, 1, 1' and 1" and one three- 
dimensional representation 3. The one-dimensional representations are given by: 

1 S = 1 T = 1 

V S = l T = e t4n / 3 = iu 2 

I" S = l T = e t2lT / 3 = cu 



The three-dimensional representation, in a basis where the generator T is diagonal, is given 
by: 



T= u> 2 , 5 = - 2 -1 2 . (81) 





The multiplication rule for triplet representations is the following: 

3 x 3 = 1 + 1' + 1" + 3 5 + 3 A (82) 

If we denote by 

a=(ai,a 2 ,a 3 ) , b=(b 1 ,b 2 ,b 3 ) (83) 

two triplets, the singlets contained in their product are given by 

1 = ( a b) = (ai&i + a 2 b 3 + a 3 b 2 ) 

V = (ab)' = (a 3 6 3 + a x b 2 + a 2 b x ) (84) 
1" = (ab)" = (a 2 b 2 + ai6 3 + a 3 6i) 

The two triplets can be separated into a symmetric and an antisymmetric part: 

3 S = (ab)s = ^(2ai6i - a 2 b 3 - a 3 b 2 , 2a 3 b 3 - a x b 2 - a 2 b±, 2a 2 b 2 - a x b 3 - a 3 bi) 

3 A = (ab) A = -(a 2 b 3 - a 3 6 2 ,ai6 2 - a 2 h,a 3 b 1 - a x b 3 ) (85) 

Moreover, if c, d and c" are singlets transforming as 1, 1' and 1", and a = (a x ,a 2 ,a 3 ) is 
a triplet, then the products ac, ad and ad' are triplets explicitly given by (aiC, a 2 c, a 3 c) , 
(a 3 d , o^c', q 2 c' ) and ( c"), respectively. Note that due to the choice of complex 

representation matrices for the real representation 3 the conjugate a* of a ~ 3 does not 
transform as 3, but rather (a*, ag, a^) transforms as triplet under A 4 . The reason for this 
is that T* = U T TU and S* = U T SU = S where U is the matrix which exchanges the 2nd 
and 3rd row and column. 



B Canonical Normalization of K and K c & Diagonal- 
ization of y\ 



We perform the following transformations on the fields present in the Lagrangian: at first 
K and K c are brought into their canonical form and then j/j is diagonalized. 

To diagonalize the hermitian matrices K and K c we apply the unitary transformations 
W and W c : 

W^KW = diag and (W C )^K C W C = diag . (86) 

Normalizing K and requires a rescaling of the fields via the real (diagonal) matrices R 
and R c : 

RW ] KWR = ± and R C (W C ) ] K C W C R C = 1 . (87) 

The fields I and e c present in the original Lagrangian (here e c stands for the three SU(2) 
singlet leptons) are expressed as: 

l = WRl' and e c = W c R c (e c ) f (88) 

so that IK I = V[RW^KWR)l' = I'll' and e c K c e c = {e c )'[R c {W c Y K c W c R c ]{e c )' = 
(e c )'l(e c y holds, yi and M. given in the basis V and (e c )' take the form: 



^/ = (e c )' T i? c (^ c ) T ?/ z Wi?r , (89) 



and 

(e c ) T Ml = (e c )' T R c (W c ) T MWRl' . (90) 

We diagonalize the resulting Yukawa couplings of the charged leptons, R c (W c ) T yiWR, by 
the usual bi-unitary transformation: 

U T [R c (W c ) T yi WR)V = ^dmg(m e ,m^m T ) (91) 

and arrive at the mass eigenbasis I" and (e c )": 

(e c )' = U(e c )" and Z' = W . (92) 

Finally, the matrix M. for the dipole moments is given as: 

(e c ) T M I = (e c )" T [U T R c (W c ) T MWRV] I" = (e c )" T M I" . (93) 

Using M. we can immediately read off the size of MDMs, EDMs and LFV processes. 

As we assume for the actual calculation of M. that all couplings involved are real, the 
matrices W, W c , U and V turn out to be orthogonal instead of unitary. Furthermore, 
we express the small parameter u in terms of t as u = xt with \x\ < 1 according to eq . 
( |T5i) and eq. ( |T8l) . We then can do the calculation in just one expansion parameter t o. 
In the course of the calculation we pose the following requirements: the kinetic terms are 
canonically normalized up to and including 0(t 5 ), yi is diagonal also up to and including 
0(t 5 ), and the matrices W, W c , U and V are orthogonal up to the same order. The 
calculations have been performed with two independent methods. 
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The different factors of t and u can be recovered in the final result by replacing x with ~ . 
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